Determination of exchange energies in the sawtooth spin ring {M075V20} by ESR 
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The magnetism of the polyoxometalate cluster {M075V20}, containing a sawtooth ring of 10 
corner-sharing triangles located on the equator of the barrel-shaped molecule, has remained de- 
batable since it is masked by contributions from impurities as well as temperature-independent 
paramagnetism. In this article we demonstrate the usefulness of ESR measurements since the 
temperature dependence of the ESR intensity can discriminate between impurity and molecular 
contributions. We determine the exchange parameters and therefore also the low-lying spectrum of 
{M075V20}, especially the low-lying singlet states which so far have been probed solely by specific 
heat measurements. 

PACS numbers: 75.10.Jm,75.50.Xx,75.40.Mg,75.50.Ee 

Keywords: Heisenberg model, Frustrated spin system. Numerically exact energy spectrum 
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I. INTRODUCTION 

The series of Keplerate molecules {Mo72Fe3o}, 
{Mo72Cr3o}, {M072V30} and {W72V30} is one of the 
beautiful creations made possible by recent developments 
of modern chemistry»i"— In these nanosized molybdate 
or tungstate-based molecules, 30 magnetic ions are lo- 
cated on the 30 vertices of an icosidodecahedron and 
are antiferromagnetically coupled, resulting in triangu- 
lar and pentagonal networks. As a result, this leads to a 
strongly frustrated cluster with a huge number of quan- 
tum states. ^^ Since the number of states, (2s -I- 1)'^", can 
be varied by substituting the magnetic ions (2^° and 6'^'^ 
for N {s — 1/2) and Fe (s — 5/2), respectively), these are 
ideal systems for studying the transitions from quantum 
to classical behavior. 

Among synthesized derivatives of the Keplerate 
clusters, the compound {M075V20} is analogous to 
{M072V30} except that the 10 V''+ ions (s = 1/2) located 
adjacent to the north and south poles of the icosidodeca- 
hedron are substituted by nonmagnetic ions.-^ Therefore, 
this compound is equivalent to a sawtooth chain of 10 tri- 
angles with a periodic boundary condition as schemati- 
cally shown in the inset of Fig. [TJ Basically, {M075V20} 
is classified as a sawtooth spin ring rather than a Ke- 
plerate cluster since it has only a partial substructure 
of {M072V30}. However, determining the exchange cou- 
plings of {M075V20} can be a good reference for the 
charactarization of its analogous compound {M072V30}. 
Moreover, sawtooth chains also belong to the class of 
frustrated antiferromagnetic spin systems with poten- 
tially very unusual magnetization curveS j^^d^ as partly 
realized in the recently investigated magnetic material 
azurite»i^ii^ 

For {M075V20}, it turns out that the product, xT, 
of the magnetic susceptibility and temperature decreases 
steadily with decreasing temperature due to antiferro- 
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FIG. 1: Temperature dependence of xr from Ref.y. The inset 
shows a scheme of the metal skeleton of {M075V20} highlight- 
ing the nearest-neighbor exchange pattern between the V*''' 
ions. 



magnetic couplings, and it was assumed in Kef's that the 
ground state is a singlet. Additionally, those authors esti- 
mated the exchange couplings between the magnetic ions 
from a theoretical fit to a reduced model that contains 
only a ring of six triangles. They concluded that their ex- 
perimental results are well reproduced when Ji = 288 K, 
J2 = 0.55 • Ji, compare inset of Fig. [H together with an 
additional coupling J3 = 0.20 • Ji between the tips of 
every second triangle.— 

Takemura and Fukumoto refined those parameters by 
using a finite-temperature Lanczos method for the full 
system of 10 coupled triangles.— They found similar ex- 
change coupling parameters, which are Ji = 388 K, 



J2 = 0.42 • Ji, and J3 = 0.21 • Ji. However, these values 
are still under discussion for various reasons. The main 
concern is the unrealistically large exchange J3 which 
should act across a distance of 11.7 A, mediated by multi- 
center exchange pathways. Another problem is given by 
the fact that the susceptibility data are superimposed by 
an unknown amount of free vanadium ions (impurities) 
as well as by temperature-independent (i.e. van Vleck) 
paramagnetism typical of polyoxometalates. 

We propose in this article that ESR is a unique method 
that overcomes these difficulties. It can separate the sig- 
nals stemming from the intrinsic {M075V20} and from 
impurities since each ESR linewidth is qualitatively dif- 
ferent. This enables us to obtain the pure magnetic re- 
sponse from {M075V20} in contrast to the susceptibil- 
ity measurements in which the intrinsic and extrinsic re- 
sponses are mixed. Since the temperature dependence of 
the ESR intensity is related to the energies of the excited 
states, we are able to obtain information about the ex- 
change couplings. The ratio of Ji and J2 determines the 
frustration of the sawtooth chain of which the density 
of low- lying singlet states is a fingerprint J^ii^ Usually 
the density of singlet-levels is deduced from specific heat 
measurements, but we will show below that ESR is also 
capable of determining the density of low-lying singlet 
states. 

The article is organized as follows. In Section [II] we 
briefly explain the experimental method. Section |llT] con- 
tains our experimental results that are compared to the 
calculated ESR intensities assuming various sets of model 
parameters. The article closes with a short summary. 



II. EXPERIMENTAL 

The Terahertz Electron Spin Resonance Apparatus in 
the Institute for Materials Research (TESRA-IMR) of 
Tohoku University has been used for the high field ESR 
measurements i^ A simple transmission method with 
Faraday configuration has been employed. We used con- 
ventional Gunn oscillators for the millimeter wave radi- 
ation, and an InSb detector for transmission detection. 
In addition, a pulsed magnetic field up to 8.5 T can be 
generated from a 90 kJ capacitor bank. Two types of 
cryostats, a conventional ^He bath type cryostat and gas- 
flow type cryostat, were used for the measurements for 
the low and high temperature ranges, respectively. Pow- 
der samples were used in this study. 



III. RESULTS AND DISCUSSION 

A. High field ESR results 

Figure [2] shows the typical temperature dependence 
of the ESR spectra for {M075V20}. The employed fre- 
quency is 190 GHz, and the temperature is varied from 
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FIG. 2: Temperature dependence of ESR spectra for 190 GHz. 
The very narrow peak observed around 6.8 T is a field marker. 
The inset shows the spectrum for 190 K. The solid and broken 
lines are fitting curves for resonance a and P, respectively (see 
text for details). 



1.5 to 250 K. As shown in the inset of Fig. [51 two absorp- 
tion lines, a broad and a sharp one, are clearly observed 
at 190 K. Hereafter, we denote the observed resonances 
as a and /? for the former and the latter, respectively. 
The tiny absorption observed at around 6.8 T is from 
DPPH, which stands for (2,2-diphenyl-l-picrylhydrazyl), 
a field marker. The intensity of a gradually increases by 
decreasing the temperature, but quickly diminishes be- 
low 80 K, which is a typical ESR behavior of the excited 
states. On the other hand, the intensity of f3 is inversely 
proportional to the temperature, following Curie's law. 
Therefore, a can be assigned to the resonance originat- 
ing from {M075V20}, and /3 from the impurities at the 
cation sites4 Since the singlet ground state is ESR silent, 
the absence of intrinsic ESR signal from {M075V20} at 
low temperature and the ESR observation of the excited 
states for relatively high temperatures suggests that the 
ground state of {M075V20} is a singlet, which is consis- 
tent with magnetic susceptibility results^"^ 

In the inset of Fig. [2] the absorption lines of a and /3 
are fitted by Gaussian curves (solid and broken curves, 
respectively), and integrated intensities for each temper- 
ature are obtained. The integrated intensity of impurities 
(i.e. of resonance /3), /imp, versus temperature is shown 
as solid circles in Fig. [3l According to the Boltzmann 
distribution, the ESR intensity of the s = 1/2 impurities 
can be written as 
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w is a coefficient which is proportional to the power 
and frequency of the radiation, A'imp is the number 
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FIG. 3: Solid circles are the integrated intensities for reso- 
nance P coming from paramagnetic impurities; they are fit- 
ted using equation ([l| (thin solid curve). Open circles are 
the normalized integrated intensities for resonance a coming 
from the intrinsic {Mo75V2o}. Thick solid and dashed curves 
are calculated curves of the normalized integrated intensity 
for various values of Ji. The calculated curves are obtained 
from the of singlet and triplet energy levels (inset) assuming 
an equilateral triangle model. The energy is normalized by 
Ji. 

of impuritieSf^ and /3 = l/ksT the inverse temper- 
ature, g is the spectroscopic splitting factor and /is 
the Bohr magneton, respectively. By using equation ([T]) 
wNi^p = 83.9 ± 7.4 is obtained from the fitting curve in 
Fig. |3] (thin solid curve). 

On the other hand, the intrinsic ESR intensity of 
the excited triplet states, /jnt, which are the transitions 
\M = -1) -^ \M ^ 0) and \M == 0) ^ \M = +1), is pro- 
portional to the radiation coefficient w and the difference 
of the population between the transition levels, i.e. 
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The factor of 2 corresponds to the square of the transition 
matrix element. ^^ When the intrinsic intensity /jnt (i.e. 
intensity of a) is normalized by wNi^p, which was al- 
ready obtained above, the unknown radiation coefficient 
w can be eliminated, and then we obtain 
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where Z{T^ B) is the partition function 
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Nint is the number of {M075V20}, and Ek and Ei are 
the energies of the singlet and triplet states, respectively. 
The normalized values of J^^^ ^^^ presented as open cir- 
cles in Fig. [21 The error bars are obtained from the 
uncertainty of the fitting curves. Using ([3]) /^„j can be 
evaluated exactly using the energy levels of the singlet 
and triplet states which can be obtained from the model 
Hamiltonian by diagonalization. This way the exchange 
constants can be determined by comparing experimental 
and theoretical values of I'^^^ . 

It is important to note that the whole procedure works 
with only singlet and triplet levels since the exchange 
constants are large and therefore states with S > 2 are 
irrelevant. 



B. Theoretical Models and analysis 

The magnetism of {M075V20} is modeled by a Heisen- 
berg Hamiltonian augmented with a Zeeman term 
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Here J^ is the exchange interaction between spins at sites 
i and j. The eigenvalues of this Hamiltonian can either 
be determined by complete matrix diagonalization or for 
the low-lying levels by the Lanczos method, ^^"^^ For our 
simulation we considered the coupling Ji on the equator 
and J2 to the tips of the triangles as schematically shown 
in the inset of Fig. [T] The low-energy part of the singlet 
and triplet sectors for J2 = Ji (i-e. equilateral triangle) 
is shown in the inset of Fig. |31 Note that energies are 
given as multiples of Ji. The lowest singlet level is de- 
generate due the frustration of the triangular magnetic 
structure.— 

From the energy eigenvalues we can calculate I^^^ . The 
ratio Mnt/Mmp — 1/2 is used since impurities of about 
2/iB per molecular unit are observed in the magnetiza- 
tion measurements. It is also important to note that we 
have taken into account all singlet and triplet levels up 
to lOEg, where Eg is the gap between the ground state 
and the lowest triplet level (see inset of Fig. |3]). In view 
of the experimental temperature range this is more than 
sufficient since the exchange parameters turn out to be 
rather large. 

The calculated //^^ curves for various Ji — J2 are 
shown in Fig. [S] The position of the maximum of the 
calculated curves changes by varying the exchange cou- 
pling constant. The same holds true for the initial rise of 
the curve which is related to the energy gap between the 
singlet and triplet states, i.e. Eg. On the other hand, the 
high-temperature behavior of J/„^ is in part related to the 
number of low-lying singlet states?^ therefore, the high- 
temperature tail is almost independent of the exchange 
parameters. For Ji = J2 = 400 K (thick solid curve in 
Fig- El), the calculated curve fits well with the initial ris- 
ing part and the tail part. However, a large difference is 



seen on the peak position. This suggests that an equi- 
lateral triangle model (J2 = Ji) is not suited for this 
system, and that an isosceles triangle model (J2 7^ Ji) 
should be considered. 
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FIG. 4: Normalized integrated intensities for resonance a 
coming from the intrinsic {M075V20} for various ratios J2/J1. 
The respective exchange parameters are shown in the legend. 

Hence, we have investigated chains of isosceles trian- 
gles with various ratios J2I J\- Figure |4] shows the curves 
for J1I J\ = 0.8, 1.2, 2.0. The absolute values of the data 
are determined by the low-temperature behavior, i.e. the 
singlet-triplet gap. It is obvious that a much improved 
fit to our experimental data is provided by assuming the 
isosceles triangle model. The agreement in the observ- 
able /(jjj appears to be better if J2 > J\. We emphasize 
that there is no need to introduce the additional exchange 
coupling J3 between tips of neighboring triangles. This 
exchange is very unlikely anyway due to the very long 
exchange pathway. 

In order to rationalize the differences between J\ and 
J2 it is helpful to look at the crystal structure of an in- 
dividual triangle of {M075V20} as schematically shown 
in Fig. [5j While the bonding length between vana- 
dium atoms along J\ and J2 bonds does not vary much, 
the V-O-Mo bonding angles do. Moreover, since the 
V-O-Mo bonding angles are larger along J2 bonds the 
Goodenough-Kanamori ruleSS suggests that the exchange 
should be stronger along this pathway, which is in ac- 
cord with our observations. Among the two good fits, 
J2/J1 = 1.2, Ji = 350 K and ,hl Ji = 2.0, Ji = 175 K, we 
tend to favor the second parameter set since this repro- 
duces the low-temperature behavior much better. The 
absolute numbers are in good agreement with other poly- 
oxometalates containing V^"'' spin centersi^i^ 

Our results suggest that the type of frustration which is 
present in {M075V20} is different from that in the original 
Keplerate molecules which are akin to the kagome lattice. 
Not only is {M075V20} a quasi one-dimensional object. 
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FIG. 5: Bond lengths and angles for an individual triangle of 

{M075V20}. 



in addition the magnetic centers are not equivalent. The 
special ratio of J2/ Ji — 2.0 relates it strongly to sawtooth 
chains with flat bands of one-magnon energies^i^^ Such 
chains are characterized by giant magnetization steps of 
50 % of the saturation magnetization. Unfortunately, the 
exchange interactions present in {M075V20} are much too 
large in order to observe such magnetization steps experi- 
mentally, even for the first one, let alone subsequent steps 
at larger fields. 



IV. SUMMARY 



In summary, we have performed ESR measurements of 
the polyoxometalate cluster {Mo75V2o}- We succeeded in 
separating the intrinsic ESR signal of the molecules and 
the signal of the magnetic impurities. Via the tempera- 
ture dependence of the integrated intensity we could de- 
duce the parameters of the underlying Heisenberg Hamil- 
tonian. We also found that the additional exchange cou- 
pling J3, which has been taken into account in previous 
studies^iii, is not necessary to explain the magnetic prop- 
erties of {M075V2o}- 

Finally, we would like to stress that ESR is a powerful 
tool to study the low energy spectrum especially of frus- 
trated magnetic systems. Since the sensitivity of ESR is 
very high, this method can in certain cases be far more 
effective than the specific heat measurements in charac- 
terizing the low-lying density of states. This study shows 
only ESR results for a single frequency of 190 GHz. It 
is of course possible to obtain more detailed information 
on the spectrum by combining results using multiple fre- 
quencies. 
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